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ABSTRACT 


As part of the Innovative Coastal-Ocean Observing 
Network (ICON), a receiver located on Sur Ridge monitored 
transmissions of low frequency tomography signals from a 
sound source on Davidson Seamount. The received signals 
were transmitted via underwater cable to the Point Sur 
Ocean Acoustics Observatory (OAO) from July 1998 through 
December 1999. Processed signals revealed a stable, 
resolvable arrival pattern. Subsequent analysis included 
forward acoustic modeling to calculate predicted raypaths. 
Observed arrivals were then associated with modeled 
raypaths, extracting observed travel times over the 17- 
month time series. Using a stochastic inverse approach, 
the extracted travel times were inverted for spatial and 
temporal variations of sound speed. Sound speed 
perturbation estimates were converted to temperature 
perturbations and compared to in situ mooring data, CTD 
transects along the acoustic path, and TOPEX/POSEIDEN 
satellite altimetry. Comparisons revealed that the 
tomographic estimate is in general agreement with the in 
situ point measurements and the altimeter data. The 
methods discussed in this paper demonstrate the application 
of ocean acoustic tomography to study temperature 
variability along the central California coast. 
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I. INTRODUCTION 


A. PHYSICAL OCEANOGRAPHY OF CENTRAL CALIFORNIA 

The California coast is an oceanographically complex 
region, little understood and inadequately sampled. 
Increased understanding of the features and dynamics of 
this region can aid fisheries and wildlife management, 
prediction and abatement of pollution and toxic 
phytoplankton blooms, atmospheric and climate change 
forecasts, and shipping and military operations (Miller, 
1999) 

Large-scale atmospheric forcing in the eastern Pacific 
Ocean consists of the North Pacific High, the Aleutian Low, 
and in summer the thermal low over the western United 
States. The North Pacific High is most intense during the 
summer months while the Aleutian Low is most intense during 
the winter months. During summer months, interaction 
between the North Pacific High and the thermal trough 
inland helps to strengthen and expand the northerly surface 
winds along the west coast of the California. This forcing 
contributes to the creation of the anticyclonic North 


Pacific Subtropical Gyre. 

The 

California 

Current System 

(CCS) forms the 

eastern 

limb 

of 

this 

gyre. 

flowing 

equatorward from 

Washington 

to 

Baja 

California. 

The CCS 

has traditionally 

been divided 

into 

three 

large- 

■scale (> 


500 km) alongshore currents: the California Current (CC) , 

the Davidson Current (DC), and the California Undercurrent 
(CUC) (Hickey, 1998). 

The CC is a surface (0-300 m deep) current, which 
carries colder, fresher subarctic water equatorward 
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throughout the year with average speeds generally less than 
25 cm s~^ (Reid and Schwartzlose, 1962) . Thus, the CC is 
characterized by a low salinity, low temperature core which 
usually lies between 300-400 km offshore (Lynn and Simpson, 
1987) . The cue flows poleward over the continental slope 
from Baja to at least 50° N with a relatively narrow width 
between 10-40 km (Hickey, 1998). The CUC has its origin in 
the eastern equatorial Pacific, and is identified by its 
warm, saline, oxygen and nutrient-poor signature. Hickey 
(1979) concluded that the location, strength and core depth 
show considerable seasonal variability and can be related 
to the seasonal variability in wind stress and curl of the 
wind stress. Peak speeds of the undercurrent are about SO¬ 
HO cm s~^, being stronger at depths of 100-300 m, and can be 
continuous over distances of more than 400 km along the 
slope (Collins et al., 1996). The DC is a seasonal 

current, flowing poleward at the surface during the fall 
and winter over the shelf from Point Conception to 
Vancouver Island. Measurements in the region have shown 
that the seasonal cycle over the slope is highly variable 
with the poleward flow maximum usually occurring in May 
(Collins et al., 1996). The reversal of winds from 
northwesterly in summer to southeasterly in winter, which 
causes downwelling at the coast, seems to be the forcing 
mechanism of this poleward surface current (Huyer et al., 
1989) . It has been suggested that the DC is a result of 
the "surfacing" of the CUC during late fall (Pavlova, 1966; 
Huyer and Smith, 1974). 

The Central California Coast, particularly in the 
vicinity of Monterey Bay, has the added complexity of 
highly variable bathymetry, which influences currents that 
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make up the CCS. This oceanographically complex region, 
which consists of variable currents, mesoscale eddies and 
upwelling events, exhibits significant variability from 
synoptic to interannual scales. It became apparent that a 
long-term, near real-time coastal observing network was 
needed to collect information on critical ocean parameters, 
for input into predictive models for the purpose of 
forecasting coastal ocean conditions. 

B. THE INNOVATIVE COASTAL-OCEAN OBSERVING NETWORK 

Funded by the National Ocean Partnership Program 
(NOPP), an alliance was formed between government, 
academic, and industrial entities to implement an 
Innovative Coastal-Ocean Observing Network (ICON) to study 
the oceanography in the Central California Coastal waters 
surrounding Monterey Bay (Paduan et al., 1999). The 
partnership involves a consortium of scientists and 
engineers from eight organizations, including the Naval 
Postgraduate School (NPS), Monterey Bay Aquarium Research 
Institute (MBARI), California State University at Monterey 
Bay (CSUMB), University of Southern Mississippi (USM), 
University of Michigan (UM), HOBI Labs, CODAR Ocean Sensors 
Ltd. (COS) , and the Naval Research Laboratory (NRL) . The 
goals of ICON were to bring together modern measurement 
technologies, to develop new technologies, and to integrate 
them within a data-assimilating coastal ocean circulation 
model. 

The principal components of the observing network 
include 1) surface current maps from shore-based high 
frequency (HF) radar installations, 2) subsurface currents. 
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temperature, salinity and bio-optical properties plus 
surface meteorological properties from several deep-sea 
moorings, 3) sea surface temperature and color from 
satellites, and 4) along-track temperature and temperature 
variances from two acoustic tomography slices through the 
region. 

The acoustic tomography component consists of two 
autonomous sound sources, one placed on Davidson Seamount 
and the other on Pioneer Seamount. These seamounts allowed 
for the placement of transmitters near the sound channel 
axis. In addition, the transmitter would be in close 
proximity to the bottom, which keeps mooring motion within 
tolerable limuts. A receiver located on Sur Ridge 
monitored the acoustic transmissions from these two sources 
continuously. The data were transmitted real time via 
underwater cable to the NPS Ocean Acoustic Observatory 
(OAO) facility at Point Sur. 

C. OCEAN ACOUSTIC TOMOGRAPHY 

Ocean acoustic tomography is a technique for observing 
the dynamic behavior of ocean processes by measuring the 
changes in travel time of acoustic signals transmitted over 
a number of ocean paths (Spindel, 1986). The word 
tomography is derived from the two Greek roots meaning "to 
slice" and "to look at." Tomography is an imaging 
technique that inverts propagation measurements through 
many sections of a volume to determine the physical 
characteristics of the interior of the volume (Medwin and 
Clay, 1998). Analogous techniques are employed in medical 
tomography (e.g.. X-rays in Computer Assisted Tomography 
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(CAT) scans) or geophysical tomography, which uses manmade 


shock waves to explore the earth's interior. Munk and 
Wunsch (1979) first proposed application of these 
techniques to the ocean as a means for monitoring ocean 
basins for mesoscale fluctuations. Ocean acoustic 

tomography uses sound energy to "look at" a "slice" of the 
ocean by measuring the travel times of various signals 
propagating along different acoustical paths through the 
volume. As acoustic energy travels along its path, sound 
speed fluctuations affect the travel time. Since sound 
speed in the ocean is a function of temperature, salinity 
and pressure, information is acquired about these 
parameters as well as current speeds and direction. These 
data are in the form of sound pulse travel time changes. 
Using these measurements of travel time perturbations, an 
estimate of the ocean structure can be constructed using 
mathematical inverse techniques. 

Ocean acoustic tomography has several advantages over 
more traditional oceanographic study methods (Chiu et al., 
1987) . An acoustic tomography system can be installed in 
the ocean as a semi-permanent, continuous, weather- 
independent observing system. The low spatial attenuation 
rate of sound and high temporal resolution allows the 
system to monitor large volumes of the ocean interior, 
sampling many different levels simultaneously. Since 
relatively few acoustic moorings are required for ocean 
acoustic tomography, costs associated with such a system 
are significantly lower than conventional "spot" mooring 
systems. Furthermore, with each additional conventional 
mooring only one piece of information is added such that a 
1:1 mooring increase to information-gain ratio occurs. In 
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contrast, the addition of one tomographic mooring adds many 
new distinct ray paths, each of which adds a piece of 
information to the system (Munk and Wunsch, 1979). 

The applicability of acoustic tomography for 
monitoring ocean variability is dependent upon the 
following four important issues: 1) stability, 2) 
resolvability, 3) identiflability, and 4) signal-to-noise 
ratio (SNR). Stability addresses the property of whether 
or not the same individual arrival exists over successive 
transmissions. Stability requires transmitter-receiver 
paths to be insensitive to changes and not fade away or 
disappear. These transmitter-receiver paths are called 
eigenrays. Resolvability requires that the arrival time 
separation between eigenrays be sufficient to resolve 
individual rays. The third issue, identiflability, 
requires that measured arrival times of eigenrays 
correspond to modeled arrival times, thus associating 
arrivals to raypaths. Lastly, adequate signal-to-noise 
ratios are required to ensure travel time precision, and 
that the signals are strong enough to be received over 
background noise. 

An ocean acoustic tomography experiment can be divided 
into two separate and distinct parts. The first is the 
forward problem, which establishes the physical 
relationship between the data and the unknown structure 
given the characteristics of the sound channel and sensor 
configuration. The second is the inverse problem, which 
demands calculation of the unknown ocean perturbations 
given raypath geometries and travel time measurements. 
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D. 


THESIS OBJECTIVES AND APPROACH 


The principal objective of this thesis is to study the 
temperature variability along the California coast between 
Davidson Seamount and Sur Ridge during the period from 30 
July 1998 to 31 December 1999. During this timeframe, the 
sound source at Davidson Seamount emitted tomographic 
signals continuously, with the exception of intermittent 
data gaps. Although the data gaps preclude the analysis of 
fast temperature fluctuations, the overall time series 
allows for an adequate analysis of mesoscale variability. 
In conjunction with mooring temperature data in the 
vicinity of the acoustic path, the Davidson Seamount to 
Point Sur transmission has allowed for an in-depth 
evaluation of the feasibility of using acoustic tomography 
in monitoring the coastal environment for mesoscale 
fluctuations. 

The analysis of the tomographic data has been divided 
into three distinct areas: 

1. Forward Modeling. Acoustic propagation modeling 
was first conducted to examine the expected ray 
arrival structure. The acoustic modeling uses a ray 
theory approach, which incorporates a reference or 
"background" sound speed profile derived from data 
acquired in July 1998 by the R.V. POINT SUR and a high 
resolution bathymetric dataset obtained from the 
National Oceanic and Atmospheric Administration 
(NOAA). To validate the model results, predicted 
arrivals were compared to observed arrivals, 
associating the observed arrivals with the predicted 
eigenrays. 
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Time Series Construction. 


Travel times for stable 


2 . 

and strong ray arrivals were extracted from the time 
series data. A 7^'^-order polynomial was used to 
interpolate through data gaps and smooth edge effects, 
while a 4*^'^-order, 10-day low-pass filter was used to 
filter out fast fluctuations such as tidal 
oscillations. 

3. Inverse Analysis. The time series of ray travel 
times were then "inverted" using a minimum mean-square 
error estimator. The temporal and horizontal 
structure of temperature along the acoustic path was 
estimated by constraining the vertical structure using 
a 3-layer model approach. Validation of the 
tomographic inverse and interpretation of the 
oceanographic variability were aided by mooring 
temperature data, CTD transects and TOPEX/POSEIDEN 
satellite altimetry. 

E. THESIS OUTLINE 

The remainder of this thesis consists of four 
additional chapters. Chapter II includes a discussion of 
acoustic tomography and oceanographic measurements as they 
pertain to this thesis. Chapter III presents and discusses 
forward modeling results, the extracted time series of ray 
travel times and ray identification. Chapter IV presents 
the vertical structure and constraint of the tomographic 
estimate, discusses development and results of the inverse 
problem, and compares the observed acoustic variability in 
relation to observed oceanographic processes. Conclusions 
are presented in Chapter V. 




II. OBSERVATIONS 


A. ACOUSTIC TOMOGRAPHY 
1. The Experiment 


Source, Receiver and Mooring Locations 



Figure 1 - Contour map depicting the orientation of 
the tomography source and receiver 

The tomography experiment extended from a transmitter 
located on Davidson Seamount northeast to a receiver placed 
on Sur Ridge. This region and location of the transmitter 
and receiver are shown in Figure 1. The transmitter and 
receiver were located at depths of 1270 and 1359 m, 
respectively. The transmitter produced phase-modulated 
signals, with a source level of 180 dB re IpPa, center 
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The 


frequency of 400 Hz, and bandwidth of 100 Hz. 
tomographic signals recorded at Sur Ridge were transmitted 
via underwater cable to the NFS Point Sur OAO. 

Transmissions were recorded from July 1998 through December 
1999. Figure 2 depicts a vertical cross-section of the 
bathymetry along the 66-km acoustic path. The bathymetry 
was interpolated from a high-resolution digital dataset 
extracted using GEODAS (GEOphysical DAta System), an 
interactive database management system obtained from the 
NOAA National Geophysical Data Center (NGDC). 


Davidson Seamount to Sur Ridge Depth Profile 



Figure 2 - Davidson Seamount to Sur Ridge Bathymetry 

In a bottom-limited configuration, ray propagation 
becomes more complex as rays interact with the bottom. 
Complicated multipath arrivals or "micro-multipaths" may 
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reach the receiver simultaneously, which exacerbates the 
problem of identifying individual observed arrivals with 
modeled ray paths. Figure 3 offers a three-dimensional 
view of the source-receiver configuration with associated 
bathymetry and example eigenray paths. 


Source/Receiver Locations, Bathymetry and Example Eigenrays 



Latitude (deg) 


Figure 3 - Sample eigenrays, Davidson Seamount source, 
and Sur Ridge receiver 

During the course of data acquisition, numerous 
outages at the NFS Point Sur OAO resulted in data gaps in 
the travel time series. Significant outages (> 2 days) are 
summarized in Table 1. 
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Data Outage 
(year/day) 

Length 

(hrs) 

Reason for Outage 

1998 226-231 

119.3 

loss of building power 

1998 318-321 

49.1 

disk storage 100% full 

1999 012-147 

3239.1 

time synchronization lost 

1999 147-167 

498 

time synchronization errors 

1999 251-297 

1104.55 

power transformer damaged 

1999 338-344 

159.48 

power outage - blown breaker 

2000 001 


all transmissions stop 

Length 

of all outages: 5356.04 hrs 

Percent of study period: 43.06% 


Table 1 - Significant outages affecting data continuity 


2. Signal Characteristics 

The traditional approach to signal design and 
processing in ocean acoustic tomography is that described 
by Birdsall and Metzger, 1986. Their method exploits phase 
modulation of the periodic signals using maximal length 
(binary) sequences (m-sequences) to achieve optimal pulse 
compression with no side lobes. Modulating the signal's 
phase with an m-sequence and processing with a correlation- 
matched filter results in a deterministic narrow pulse, 
which leads to a signal processing gain of 101og(511) 
(where 511 is the m-sequence length). Transmitting this 
sequence repeatedly leads to additional signal processing 
gain of 101og(35) at the receiver by coherently averaging 
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the consecutive sequences (35 m-sequences transmitted). 
Source signal characteristics are summarized in Table 2. 


Source Characteristics - Davidson Seamount 

Source level 

180 dB re IpPA at 1 m 

Center frequency 

400 Hz 

Full bandwidth 

100 Hz 

m-sequence law (octal) 

1021 

m-sequence length 

511 digits 

(signal processing gain of 27 dB) 

Digital width (resolution) 

4 cycles = 10 msec 

m-sequence period 

5.11 sec 

m-sequences transmitted 

35 sequences 

(signal processing gain of 15 dB) 

Transmission length 

178.85 sec 


Table 2 - Characteristics of the 400 Hz tomography signal 


The sound source was programmed to transmit signals at 
two different rates, for the purpose of resolving 
oceanographic variability. From 30 July to 28 August 1998 
the source transmitted every 30 minutes, the objective 
being to resolve high-frequency oscillations such as 
internal tides. From 28 August 1998 thereafter, the source 
transmitted every 12 hours to resolve longer-term 
variations. At the Point Sur OAO, two different sampling 
rates were utilized. A sampling rate of 1000 Hz was 
initially utilized. Even though this sampling rate met the 
Nyquest criterion, it proved insufficient to produce an 
improved SNR through coherent averaging of multiple arrival 
structures. The sampling rate was changed to 2000 Hz on 29 
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August 1998 and continued at this rate through December 
1999. 

One of the fundamental issues in ocean acoustic 
tomography is the ability to record time accurately. A 
SeaScan temperature-compensated quartz crystal clock 
provided timing control at the source. The receiver used a 
GPS trigger to start data acquisition, with a crystal 
oscillator determining the sample rate. Significant rate 
variability was observed after power outages. This led to 
the calculation of drift rates, which are summarized in 
Table 3 along with additional timing inconsistencies 
determined during the July 1998 to December 1999 time 
series . 


Timing Error 

Data Set 

Correction 

Clock offset 

Entire time series 

-0.098307 sec 

File-by-file drift 

1998-1000 Hz 

+0.4548 sec/day 

File-by-file drift 

1998-2000 Hz (242-317) 

+0.4619 sec/day 

File-by-file drift 

1998-2000 Hz (318-365) 

+0.3629 sec/day 

File-by-file drift 

1999-2000 Hz 

+0.4012 sec/day 

Source drift 

Entire time series 

+0.00126 sec/day 

Leap Second 

1999-2000 Hz 

+1.0 sec 


Table 3 - Timing Corrections 


3. Arrival Structure 

A correlation-matched filter was applied to the data, 
with the output corresponding to a sum of pulses arriving 
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from the multiple eigenray paths. An improved signal-to- 
noise ratio was obtained through coherent averaging. Each 
transmission was processed and inspected individually. 
Figure 4 depicts a typical arrival pattern of one 
transmission while Figure 5 shows an example waterfall plot 
of transmissions from 30 July 1998 to 12 January 1999. The 

waterfall plot is useful for examination of arrival 

stability, and assists in determination of whether the 
arrival peaks are resolvable over the entire time series. 

Once all transmissions were processed, the travel 
times of the arrival peaks were extracted. Travel times 

within a window were picked for peaks that exceeded a 
median amplitude threshold value for that window. Since 
the time series covered over 17 months, the travel time of 
an individual arrival can change more than the separation 
between individual arrivals. This led to using multiple 
windows and adjusting the window accordingly to keep the 
arrival time peaks within the limits of the window. 

Applying a 7th-order polynomial to fill in small data gaps 
and smooth edge effects further refined the arrival 
structure. High frequency oscillations were removed using 
a 10-day low pass filter. Results are shown in Figure 6 as 
a dot plot (the two large gaps in the dot plot correspond 
to outages addressed in Table 1) . The dot plot allows 
identification of individual rays at the receiver location 
and match these to the modeled arrivals, which will be 
discussed in Chapter 3. 
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Normalized Anplitude 



Figure 4 - Example of coherently averaged signal 


OlDserved Arrival Structure (July 30, 1998 - January 12, 1999) 



380 


360 


220 


4S 46.5 46 
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200 


Figure 5 


Waterfall of observed arrival structure 
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Yearday - Reference 1998 























Arrival Time Dot Plot 



Figure 6 - Travel time dot plot of observed arrivals 

B. OCEANOGRAPHIC DATA 
1. Shipboard Data 

Acoustic modeling of the received signal requires 
creation of a "background" or "reference" ocean. One input 
required by the model is a reference sound speed profile. 
During July 1998, an oceanographic cruise aboard the R.V. 
POINT SUR provided CTD data at 9 stations along a path 
between Davidson Seamount and Sur Ridge. Based upon the 
CTD data, a mean sound speed profile was derived (Figure 
7) . This mean profile represents the reference ocean in 
forward modeling, and serves as an initialization for the 
inverse. 
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Figure 7 - Mean sound speed profile and reference 
sound speed field (upper 1000 m) 


2. Moored Temperature Data 

During the acoustic transmissions, oceanographic 
moorings M2 and M4, (deployed by MBARI and NFS, 
respectively) provided in-situ temperature data. The M2 
time series covered the entire period of transmissions, 
recording temperature every 10 minutes at depths of 10, 20, 

40, 60, 80, 100, 150, 200, 250 and 300 m. M4 included an 

additional temperature measurement at 350 m, but was only 
operational for the last four months of transmissions 
(September through December, 1999). Geographic locations 
of M2 and M4 are depicted in Figure 1. For an overall 
sense of ocean variability throughout the time series, M2 

temperature variations are depicted in Figure 8. 
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M2 Tenderature Data: August 1998 thru December 1999 



Figure 8 - Temperature variation at M2 


19 



























THIS PAGE INTENIONALLY LEET BLANK 


20 


III. FORWARD MODELING 


A. HAMILTONIAN RAY TRACING 

Ray theory allows the study of sound propagation 
through a medium whose refractive index can vary in space. 
Raytracing provides a graphical representation of the 
trajectories taken by sound energy through the medium. 
According to ray theory, the signal received is the sum of 
many different arrivals, each following different paths, or 
eigenrays, from the transmitter to receiver. Pressure 
disturbances at the receiver are associated with individual 
eigenrays, which have different phase shifts, time delays, 
and amplitudes. The received signal can be written as: 

r {t) = 2^ ans {t - tn)e d) 

where r(t) is the complex envelope of the received signal, 
s (t) is the complex envelope of the emitted signal, fo is 
the carrier frequency of the transmission, and tn, an, and 
On are the time delay, amplitude modification and phase 
shift of each eigenray, respectively. Consequently, 
forward modeling of the received signal demands raytracing, 
eigenray search, and determination of tn, an and On. 

The three-dimensional raytracing program HARPO 
(Hamiltonian Acoustic Raytracing Program for the Ocean) was 
used to calculate raypaths. HARPO traces each raypath by 
numerically integrating Hamiltons's equations of motion in 
spherical coordinates with a different set of initial 
conditions. In the high-frequency limit, sound waves 
behave like particles and travel along rays, according to 
the equations that govern changes of position and momentum 
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HARPO 


in a mechanical system (Jones et al., 1986). 

requires input of a model ocean by the user, which consists 
of a continuous representation of the sound speed field and 
continuous two-dimensional representation of the upper and 
lower reflecting surfaces. The upper and lower reflecting 
surfaces are the ocean surface and ocean bottom, 

respectively. Continuous treatment of sound speeds and 
bathymetry eliminates the problems of false caustics and 
discontinuous raypath properties (Jones et al., 1986). 
Since its inception, HARPO has undergone various 

modifications and improvements. Newhall et al. (1990) 

created supplementary peripheral routines for efficient 
eigenray finding and interpolation of gridded sound speed, 
current, and bathymetric data. Bathymetric and mean sound 
speed profiles used in HARPO for this study are shown in 

Figures 2 and 7, respectively. 

Within HARPO, a fan of rays from -15° to +15° was 
traced at a resolution of 0.002°. Rays greater than 15° and 
less than -15° were omitted due to numerous bottom 
interactions, which lead to a significant energy loss at 
the receiver. 

B. ARRIVAL STRUCTURE 

The output of HARPO includes eigenray geometry and 
travel time. Processing of the HARPO output was 

accomplished using the program "ray2db" developed by Chiu 
(1994). This program searches for eigenrays, computes 
travel times, calculates phase shifts as a result of bottom 
and surface reflections, and estimates signal losses due to 
raytube spreading. Program "ray2db" operates on the HARPO 
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output file utilizing input parameters set by the user in 
the initialization program, "ray2dbstart". Input 
parameters include source and receiver separation, location 
and depth, source signal characteristics, and boundary 
condition parameters. Boundary conditions affect both the 
magnitude and phase of the eigenray arrivals through 
surface and bottom reflection. A root-mean-square sea 
surface height of 1 m represented the surface boundary 
condition. At the bottom boundary, complicated 
interactions of sound waves with the upper sea floor layer 
require geoacoustic modeling to compute general values and 
restrictive parameters for various marine sediment 
properties. These sediment properties vary between the 
source and receiver. At Davidson Seamount, the bottom is 
composed of rock outcroppings interspersed with coarse 
grained sediment, while at Sur Ridge the bottom is 
primarily coarser, glauconite-rich sediment (Gabriel, 
2001) . In between, the upper layer is generally comprised 
of fine to medium grain sediment. Based upon sediment 
chart and geoacoustic parameter tables (Hamilton, 1980), 
input parameters were determined for the bottom boundary 
(Table 4) . 


Average sediment sound speed 

1650 m/s 

Average sediment density 

1600 kg/m^ 

Sediment attenuation rate 

0.04 dB/m/kHz 


Table 4 - Geoacoustic Parameters 

Given the HARPO output along with input parameters, 
"ray2db" computes the envelope of the received signal. The 
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Depth (km) 


modeled signal included 113 eigenrays, which are shown in 
Figure 9. The modeled arrival structure is shown in Figure 
10 as a stem plot of the 113 individual arrivals. 


Eigenray Tracing from Davidson Seamount to Sur Ridge 



Figure 9 - Eigenray Geometry 










Eigenray Arrivals 



Travel time in seconds 


Figure 10 - Stem plot of individual eigenray arrivals 

The original modeled arrival time included a bias 
error associated with the actual source/receiver 
separation. An iterative least squares inversion was used 
to estimate the range error. A final range error 
estimation of 0.2496 km resulted in a travel time change of 
-166 ms. This estimated travel time change was then 
applied to the modeled rays, which assisted in identifying 
observed arrivals with model results. 


25 




































c. 


RAY IDENTIFICATION 


Identification of ray arrivals was initially conducted 
by comparing the stem plot of the modeled eigenray arrivals 
(Figure 10) to the observed arrivals, waterfall plot and 
dot plot (Figures 4, 5 and 6, respectfully) . To further 

identify rays with observed arrivals, additional plots of 
the modeled arrivals were created. First, travel times 
versus launch angles were plotted, which revealed that each 
of the arrivals not only is associated with individual 
eigenrays but groups or "bands" of eigenrays with each band 
representing a set of micro-multipath that reach the 
receiver at the same time and sample the same space in the 
ocean. Next, a series of ray diagrams were constructed 
associating each eigenray band with each isolated arrival. 
Travel times versus launch angle are shown in Figure 11, 
while the geometry of the eigenray bands are shown in 
Figure 12 . 


Axrzr±v’a.l T±me versus Ija.uncli Ang^le (Modeled) 



Figure 11 - Travel time versus launch angle (modeled) 
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Arrival Bands 1-10 


(a) 



(e) 



<b) 



(f) 
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(h) 
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Figure 12 - Eigenray arrival bands 
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Using a combination of launch angle and arrival time a 
relationship was established between eigenray bands and 
observed arrivals. Each band of eigenrays was further 
examined to assist in identifying specific observed 
arrivals with modeled arrivals. Ten observed arrivals were 
associated with the 10-eigenray bands shown in Figure 12. 

A band of nine eigenrays was identified in the modeled 
structure, but could not be identified as stable and 
resolvable arrivals in the observed data (e.g., a "ghost" 
arrival). Two additional bands of eigenrays were not 
identified in the observed data (Figure 13) . Attenuation 
due to numerous surface and bottom bounces and 
deconstructive interference between the micro-multipaths 
may have contributed to a weak signal, unidentifiable in 
the presence of noise. These two unidentified rays are 
apparent in Figure 15, where their sound pressures are 
extremely low. 


Unidentified. Ajririvals 




Figure 13 - Unidentified arrival bands 
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Once all modeled eigenrays bands were associated with 
particular observed arrivals (or were deemed 
unidentifiable) , a single ray was chosen from each band to 
represent the path geometry of the band in the inversion. 
The eigenray geometry representing the 10 bands is depicted 
in Figure 14 while the modeled arrival structure is shown 
in Figure 15. 


Final Eigenray Geometry for Inverse Input 



Figure 14 - Final 10 Eigenray Paths 
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Modeled Signal Arrival Pattern 



Figure 15 - Final modeled arrivals 

D. TRAVEL TIME VARIABILITY 

Sound speed in the ocean is a function of temperature, 
pressure and salinity, with temperature dominating in the 
upper 1000 m where the gradient is strongest. As the 
temperature increases, sound speed increases thus travel 
time decreases. Likewise, as the temperature decreases 
sound speed decreases and travel time increases. Known 
ocean processes can be directly related to temperature 
changes, which can then be observed as travel time changes. 
The principle objective of this experiment is to study the 
low frequency temperature variability between Davidson 
Seamount and Sur Ridge. Significant sources of low 
frequency oscillations in the ocean include seasonal or 
inter-annual change, synoptic scale events, and mesoscale 
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variations. Seasonal oscillations can be clearly seen in 
both the mooring and travel time data, which are 
characterized by cooling of the surface waters in winter 
and warming in summer. Synoptic scale events, such as 
upwelling and downwelling, are also evident in the mooring 
and travel time data. These types of events can be 
associated with northwest winds off Central California 
resulting in the upwelling of cool, nutrient rich waters, 
and south-southeast winds which result in a downwelling of 
warmer surface waters. Although present throughout the 
year, mesoscale variations may or may not be apparent in 
the mooring and travel time data simultaneously. These 
events can be much more localized, identified at one 
location (i.e., M4) and not the next (i.e., Davidson 
Seamount). 

Semi-diurnal and diurnal variations, such as ocean 
tides and fluctuations in the mixed layer depth, were 
observed in the mooring and travel time data. These higher 
frequency oscillations were not the subject of this 
experiment and were filtered out using a order, 10-day 

low pass filter. 
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IV. TOMOGRAPHIC INVERSE 


A. DEVELOPMENT OF THE INVERSE PROBLEM 


Qualitatively, the inverse problem is relatively 
straightforward - infer the unknown ocean structure given 
the measurements of acoustic travel time along different 
raypaths. Quantitatively, the mathematical formulation 
requires some forethought. To apply inverse theory, a 
valid mathematical model must be established relating model 
parameters to measured data. A general mathematical 
expression used for ocean observing systems is the Fredholm 
integral equation of the first kind: 

d{x') = ^ g{x\x)f {^dx +e{x') (2) 

/ 


where d (x') is the data observed at some x', g(x'',x) is the 
"model" or physics of the problem, f (x) is the unknown 
function, and e (x') is the noise contamination. As opposed 
to a discrete observation, equation (2) represents a 
continuous observation over the space x' , and linearly 
relates the data to the unknown. In this experiment the 
observed data is the perturbations of the travel times of 
the resolved rays, the model is the raypath output from 
HARPO, and the unknown is sound speed change. 

Relating sound speed change to travel time change, the 
above equation can be expressed as: 


St i = 


I (-) S c(^)ds + e(^) 

path] C f X ) 


(3) 
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where St j is the change in travel time due to the sound 

speed perturbation 5c along path j, c is the reference 
sound speed field, 5c is the unknown variable which is the 
departure from the reference sound speed field, x is the 
position vector and s is the arc length along the path j 
connecting the source and receiver. Contamination by noise 
is represented by e, and includes both measurement and 
model error. 

The approach used will be a stochastic inverse method 
(Chiu et. al, 1994j, which is based on a spectral 
decomposition of the sound speed perturbations and a 
minimization of an objective function. This method treats 
the unknowns as random variables, and requires 
specification of statistical information on the signal and 
noise. Specifying ocean decorrelation scales, a noise 
variance and a solution variance the method provides an 
optimal estimate of the solution, which has minimum mean 
square errors. To aid in oceanographic interpretation, the 
method also gives solution error and resolution estimates. 
Empirical Orthogonal Function Analysis (EOFA) on 
hydrographic data is employed to investigate the vertical 
structure for constraint in the tomographic estimate. 

B. VERTICAL STRUCTURE AND CONSTRAINT 

EOFA is a common statistical method for analyzing 
oceanographic variability occurring across a spectrum of 
spatial and temporal scales. This method is used to 
redistribute the variability of a large set of variables to 
a much smaller set, which contains most of the original 

variance. In other words, EOFA estimates the spatial 
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variability and the time series coefficients modulating the 
temporal variability from the variance-covariance matrix of 
the observed data. The input is the observed data, while 
the output consists of three types of results, 1) principal 
components (PCs), 2) eigenvectors (EOFs), and 3) 
eigenvalues. The PCs, plotted as a time series, quantify 
the overall strength of the associated EOF pattern over 
time. The relative relationships between EOF points (or 
"modes") remain the same, but the absolute magnitude of the 
pattern changes with time. 

An EOFA was constructed from the M2 mooring data. 
Modes 1-3 are shown in Figure 16. The analysis revealed 
that mode 1 was dominant, accounting for 54% of the 
variance. Modes 2 and 3 contained 22% and 8% of the 
variance, respectively. Simulation studies show that they 
do not contribute significantly to travel time change. 

Additional deep CTD data analysis revealed the 
possibility of a phase change at 450-500 m and significant 
perturbations below 500 m. To account for the possible 
phase change and deep perturbations in the inverse, two 
additional layers were included. One can think of the 
intermediate (300-500 m) and deep (500-2000 m) layers as 
two additional EOFs with their unknown coefficients 
corresponding to layer-averaged perturbations. Figure 17 
shows the 3-layer model that was used to constrain the 
vertical structure in the tomographic inverse. 
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Depth (km) Depth (m) 


M2 Modes 1-3 



Figure 16 - M2 vertical structure 


Mod.e 1 a.nd. La.yer-Averaged. Modes 



Mod.a.1 Ampli tilde 


Figure 17 - Modes used for vertical constraint 
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c. 


INVERSION 


The change in sound speed is a function range, depth 
and time. This variation in sound speed can be described 
by a linear combination of n EOFs, i.e., 

n 

Sc{x,z,t) = '^ai{x,t)fi{z) (4) 

where a(x,t) is the time and range varying amplitude of the 
EOF and fi(z) is the i^'^ mode. In this description, the 
unknown becomes the coefficients of the EOFs. The 

horizontal variation of ai(x,t) was discretized with a grid 
spacing of 1 km. The vertical structure was fixed with the 
first mode and two layer-averaged modes. The change in 
travel time is referenced to the first travel time 
measured, which coincides with the July 1998 CTD transect. 
An iterative solution was employed to refine the a priori 
statistical information, specifically the rms values of 
coefficients of perturbation EOFs and noise variance. The 
inverse method processes travel times for each iteration, 
giving daily estimates of sound speed perturbation. The 
solution was deemed valid when statistics of the final 
solution and data residual were consistent with the assumed 
values for the statistical parameters. 

Several iterations of the inversion were used to 
refine the initial estimates. The input statistical 
parameters for the final iteration were 64 (ms)^ for the 

noise variance, 26 km for the horizontal decorrelation 
length, 25 m^/s^ for the solution variance of mode 1, 0.25 

rc? / for the solution variance of the intermediate depth 
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layer, and 1 m^/s^ for the solution variance of the deep 

perturbation layer. 


D. RESULTS AND DISCUSSION 

1. Resolution and Mean Square Error of the 
Tomographic Estimate 

Horizontal resolution lengths of the tomographic 
estimate are depicted in Figure 18, and represent the 
smallest ocean feature that can be seen by tomography. 
Minimum resolution lengths of about 30-35 km occur for mode 
1 and the deep perturbation at about 45 km in range, while 
for the intermediate mode the minimum resolution length of 
39 km occurs at a range of 30 km. At ranges between 0 and 
25 km, the density of raypath crossings is minimal, 
resulting in a poor resolution. Thus, in the first half of 
the transmission path, the tomographic system is unable to 
resolve ocean features that are shorter than the acoustic 
path. This poor resolution is also apparent near the 
receiver, where again the density of raypath crossings is 
minimal. 

The mean square error (in percentage) or uncertainty 
in the estimate is shown in Figure 19. The estimated mode 
1 coefficient has an average mean square error of less than 
43%. For the intermediate depth layer the uncertainty in 
its coefficient estimate is 100%, which is due to the fact 
that the corresponding travel time signal is much less than 
the travel time noise. The deep perturbation layer has an 
average uncertainty of 50%. 
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Mean Square Error (%) 


Hoirizon'ta.J. ResoJ.u't±on Iieng'ths 



Figure 18 - Horizontal resolution lengths of the 
tomographic estimate 


Eirr-or- Distribution 



Figure 19 - Error distribution in physical space 
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2. Interpretation of Tomographic Estimate 

Depth-averaged temperature perturbation estimates for 
mode 1 and the deep perturbation layer are shown in Figures 
20 and 21, respectively. Both estimates show an initial 
zero perturbation at 1998 yearday 210, corresponding to 
initialization of the inverse to the first day travel time. 

The mode 1 perturbation estimate was examined first. 
From 1998 yearday 210 to 375, the estimate showed a gradual 
cooling followed by a stronger cooling in the second half 
of the acoustic path. From 1998 yearday 575 to 735, the 
estimate showed cooling followed by slight warming, 
followed by stronger cooling. Once again the cooling is 
more defined in the second half of the acoustic path. This 
coincides with Figures 18 and 19, which show acceptable 
resolution and mean square error in this area. Although 
the error is acceptable in the first half of the acoustic 
path, the horizontal resolution indicates that the inverse 
was unable to adequately resolve ocean features smaller 
than the acoustic path. 

The deep perturbation estimate showed gradual cooling 
followed by stronger cooling up to 1998 yearday 375. From 
1998 yearday 535 to 735, the estimate revealed a strong 
warming followed by gradual and then deeper cooling. The 
error associated with the deep perturbation is marginal at 
best, with resolution lengths shorter than the acoustic 
path only between 20-55 km in range. 
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Depth Averaged Teroperature Perturbation: Mode 1 
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Figure 20 - Depth-averaged temperature perturbation 
estimate (mode 1, 0 - 300 m) 
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Depth Averaged Terrperature Perturbation: Deep Layer 
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Figure 21 - Depth averaged temperature perturbation 
estimate (deep mode, 500 - 2000 m) 
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3. Comparison to Moored Temperature Data 

Figure 22 shows M2 mode 1 coefficients along with 
range-averaged perturbation estimates for mode 1 and the 
deep perturbation layer (July 1998 - January 1999). The M2 
perturbations signify a change from the M2 mean, while the 
range average perturbations signify a change from the first 
day. It is readily apparent that M2 mode 1 and the 
tomographic solution for mode 1 are in good agreement. 
Both show an overall cooling of similar magnitude in the 
upper layer where the mode reaches a maximum. The range- 
averaged solution for the deep perturbation shows a gradual 
cooling over the time series. 
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Figure 22 - Comparison of M2 mode 1 coefficients and range 
averaged tomographic estimate (July 1998 to January 1999) 
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Figure 23 shows M2 and M4 mode 1 coefficients along 
with range-averaged perturbation estimates for mode 1 and 
the deep perturbation layer (June 1999 - December 1999). 

The M2 and M4 perturbations signify a change from the M2 
mean, while the range average perturbations signify a 
change from the first day. The overall warming and cooling 
trends of the mode 1 tomographic estimate are in general 
agreement with the mooring coefficients, except that the 
magnitude of the tomographic estimate is much larger than 
the M2 and M4 coefficients. This may be due to a localized 
cooling event along the acoustic path. The deep 

perturbation estimate shows slight warming followed by 
gradual cooling over the time series. 
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Figure 23 - Comparison of M2 and M4 coefficients, and range 
averaged tomographic estimate (June 1999 to December 1999) 
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4. Comparison to CTD Transects 

As a consistency check on the tomographic estimate, 
temperature data from CTD transects along the acoustic path 
were examined. For comparison purposes, the July 1998 
reference section is shown alongside the July 1999, August 
1999 and October 1999 sections (Figures 24, 25 and 26, 
respectively). To assist in examining the deep 
perturbation layer, deep mooring data from a current meter 
at M2 and a CTD at the source are shown in Figure 27. 

The mode 1 perturbation estimate (Figure 20) was 
examined first. Comparing the July 1998 (1998 yearday 195) 
and July 1999 (1998 yearday 560) temperature sections 
(Figure 24), significant cooling was visible throughout the 
upper 300 m. Overall, this compares well with the 
perturbation estimate. Localized events along the acoustic 
path were further examined, specifically along the second 
half of the acoustic path. During July 1999 (1998 yearday 
560), the perturbation estimate showed regions of warmer 
water surrounded by cooler water, corresponding to positive 
and negative perturbations, respectively. Upon inspection 
of the July 1999 (1998 yearday 560) temperature section, 
significant warming was evident, particularly at depths of 
200-300 m near the receiver. This warming signature may be 
due to an increase in strength of the California 
Undercurrent as the "spring transition" relaxes. The 
August 1999 (1998 yearday 590) temperature section (Figure 
25) showed a reversal of the isotherms near the receiver. 
This cooling signature can be seen as "more negative" 
perturbations in the tomographic estimate, which may be 
indicative of a relaxation of the California Undercurrent. 
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The October temperature section (Figure 26) showed further 
overall cooling across the 0-300 m layer, which corresponds 
to the stronger negative perturbations across the gap 
(October 1999 (1998 yearday 640) falls within the data 

gap) . 

The depth-averaged perturbation estimate for the deep 
layer was examined next. Comparing the July 1998 and July 
1999 temperature sections (Figure 24) between 500 and 1000 
m, a slight warming or pushing down of the isotherms was 
visible in the first 30 km. This was apparent in the 
perturbation estimate, which showed warming over the depth 
averaged layer. Analysis of deep temperature data at the 
source (Figure 27) was consistent with initial seasonal 
cooling followed by overall warming. The August 1999 (1998 

yearday 590) temperature section (Figure 25) showed 

additional warming in the deeper layers while the estimate 
revealed a gradual cooling over the depth-averaged layer. 
The October 1999 section (Figure 26) showed a reversal of 
the isotherms. A data gap in the estimate precludes 
comparison with tomography, however, this cooling was 
apparent later in the estimate and was consistent with the 
seasonal cooling visible in the deep mooring data. In the 
second half of the acoustic path, the July 1999 (1998 
yearday 560) temperature section revealed a dominant 

cooling event followed by slight warming near the receiver, 
while the August 1999 (1998 yearday 590) section revealed 

warming followed by slight cooling. The overall cooling 
signature was apparent in the perturbation estimate, 
however, the slight warming signature at the receiver was 
not. The October 1999 (1998 yearday 640) temperature 

section showed significant cooling along the second half of 
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the acoustic path, which can be interpolated in the 
estimate across the data gap. In addition, the seasonal 
cooling was readily apparent in the deep mooring data. 
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Figure 24 - Temperature data from July 1998 and July 1999 
transect along acoustic path 
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Figure 25 - Temperature data from July 1998 and August 1999 

transect along acoustic path 
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Figure 26 - Temperature data from July 1998 and October 
1999 transect along acoustic path 
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Figure 27 - Deep mooring temperature data at source and M2 

5. Comparison to TOPEX/POSEIDEN Altimetry 

TOPEX/POSEIDEN sea surface height anomalies in the 
vicinity of the acoustic path were examined over the 
tomographic time series. The altimetry data (Eigure 28) is 
in general agreement with the mode 1 perturbation estimate 
(Eigure 20) . Tomography shows a warm ocean followed by 
strong cooling, consistent with the change in sea surface 
height. An exception is near the source, where the 
altimetry shows stronger cooling. Erom 1998 yearday 550 
and later, tomography shows a cooler ocean. A noticeable 
correlation was observed in the vicinity of 1998 yearday 
675, where the estimate shows a strong cooling, which is 

consistent with the lower sea surface height shown in the 
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altimetry data. The high sea surface height visible in the 
altimetry in the vicinity of 1998 yearday 600-650 is not 
apparent in the tomographic data. This is most likely due 
to late summer surface warming, which was not detected by 
the depth-averaged tomography. 

Topex/Poseiden Sea Surface Height 
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Figure 28 - TOPEX/POSEIDEN sea surface height 
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V. CONCLUSIONS 


The objectives of this thesis were to examine 
temperature variability between Davidson Seamount and Sur 
Ridge using ocean acoustic tomography. The study included 
solving the forward problem, identifying modeled rays with 
observed arrivals, and inverting travel time perturbations 
using a stochastic inverse method. In solving the forward 


problem. 

a ray- 

■theory 

approach 

was 

used to model the 

arrival 

pattern . 

Modeled rays 

were 

associated 

with 10 

observed 

arrival 

peaks, 

and the 

time 

series of 

observed 


travel times were extracted. The vertical structure was 
constrained using a 3-layer model consisting of a dominant 
surface intensified mode, an intermediate depth layer and a 
deep perturbation layer. An iterative solution was used to 
refine the estimated statistical parameters, to arrive at a 
best estimate of the temperature perturbation field. Based 
on the estimated perturbations, coastal ocean variability 
along the transmission path was examined. 

The following summarizes major findings of this 
thesis: 

1) Horizontal resolution was poor near the source and 
receiver. Since eigenray geometry determines the 
horizontal and vertical resolving power, a lack of 
raypath crossings results in poor resolution. 
Higher resolution could be obtained by locating the 
source and receiver closer to the sound channel 
axis . 

2) The vertical structure was not well represented by 1 

mode. A 3-layer model was used to constrain the 
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vertical in the tomographic estimate. An EOFA on M2 
hydrographic data resulted in a surface dominated 
mode, which represented the upper 300 m. Additional 
deep CTD analysis revealed a possible phase change 
at 450-500 m and perturbations below 500 m. To 
account for possible phase changes and deep 
perturbations in the inverse, two additional layers 
(300-500 m and 500-2000 m) were included as two EOFs 
with their unknown coefficients corresponding to 
layer-averaged perturbations. 

3) The tomographic estimate was compared to in situ 
mooring data, CTD transects along the acoustic path, 
and TOPEX/POSEIDEN satellite altimetry. Comparisons 
show that the tomographic estimate has high enough 
spatial and temporal resolution to distinguish 
seasonal and mesoscale oscillations. The time 

series was dominated by significant cooling and 
warming events, which may have masked additional 
mesoscale fluctuations. 
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